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ABSTRACT 

We revisit the cloud-in-cloud problem for non-Gaussian density fluctuations. We show that the extended Press- 
Schechter (EPS) formalism for non-Gaussian fluctuations has a flaw in describing mass functions regardless of 
type of filtering. As an example, we consider non-Gaussian models in which density fluctuations at a point obeys 
a distribution with v degrees of freedom. We find that mass functions predicted by using an integral formula 
proposed by Jedamzik, and Yano, Nagashima & Gouda, properly taking into account correlation between objects 
at different scales, deviate from those predicted by using the EPS formalism, especially for strongly non-Gaussian 
fluctuations. Our results for the mass function at large mass scales are consistent with those by Avelino & Viana 
obtained from numerical simulations. 

Subject headings: cosmology: theory — galaxies: formation — galaxies: mass function — large-scale structure 
of universe 



L INTRODUCTION 

How many objects with mass M are there in our universe? 
This question has been one of main interests in the field of cos- 
mological structure formation. Formation process of cosmolog- 
ical objects such as galaxy clusters and galaxies is well under- 
stood qualitatively in the context of the hierarchical clustering 
scenario based on a cold dark matter (CDM) model. In order 
to compute the number density of collapsed objects or mass 
function, one must deal with gravitational non-linear growth of 
small density perturbations. The most direct way is to perform 
A^-body simulations. However, performing simulations for a 
large number of models on wide range of scales is a very diffi- 
cult task because of the limit of the computation time and avail- 
able amount of memory. Therefore, it is of great importance 
to derive analytic formulae that accurately describe the result 
of A^-body simulations. Among them, the Press-Schechter for- 
mula (Press & Schechter 1974; hereafter PS) has been a most 
successful one and applied to a wide class of structure forma- 
tion models. 

Nevertheless, the PS formalism has a flaw in describing the 
number of collapsed objects. Because the underdense region of 
smoothed density fluctuations is not taken into account in the 
formalism, integration of the mass function over whole range 
of mass does not yield a unity. Even if the density fluctuation 
smoothed on mass scale M is less than a critical density fluctu- 
ation (5c, there is a chance that the density fluctuation smoothed 
on larger mass scale M' > M is larger than 6c- This is called the 
"cloud-in-cloud" problem. Press & Schechter (1974) simply 
multiplied the mass function by the "PS fudge factor of two"in 
the case of Gaussian random fields. 

The cloud-in-cloud problem for Gaussian random fields have 
been partially solved by Peacock & Heavens (1990) and Bond 
et al (1991) using the so-called excursion set formalism and 
by Jedamzik (1995) and Yano, Nagashima & Gouda (1996) 
using an integral equation. Consider a density fluctuation 5m 
smoothed on mass scale M at a given point. Then one can 
regard a sequence of density fluctuations (5m, , 6mi ,■■ ■ in de- 
scending order of M as a trajectory of a "particle." For fluc- 
tuations smoothed by the sharp fc-space filter, each trajectory 



is described by a Markovian random walk of (5 as a function 
of "time" M. Then we analytically obtain the PS fudge fac- 
tor of two. However, for fluctuations smoothed by other filters, 
no analytic result is known, since the correlation between fluc- 
tuations with different scale renders the motion of a particle 
non-Markov process. 

In contrast, the cloud-in-cloud problem for non-Gaussian 
density fluctuations had not been explored until recently. As 
is well known, a number of theoretical "unstandard" models in- 
cluding cosmic string models, texture models, multiple fields 
models, and so on predict that the primordial fluctuations are 
not Gaussian. Because the number of collapsed dark matter ha- 
los at early epoch, e.g. clusters at z 1 or galaxies at z > 5, de- 
pends sensitively on the tails of distribution function of initial 
density fluctuations, even a small deviation from Gaussianity 
would cause a noticeable change in the statistical property of 
those high-redshift objects. To make predictions on the num- 
ber count of these rare objects, it is of crucial importance to 
investigate how the PS formalism is extended to models with 
non-Gaussian initial conditions. 

Recently, there has been some progress on this issue based 
on the so-called extended Press-Schechter (EPS) formalism in 
which the PS fudge factor is again assumed to be a constant. 
For Gaussian fluctuations smoothed by the sharp ^-space filter, 
the assumption is correct because of the nature of the Marko- 
vian random walk of 5. However, in other cases, it is not 
known whether such assumption is correct, especially for non- 
Gaussian models. In the context of formalism developed by 
Jedamzik (1995) and Yano, Nagashima & Gouda (1996), the 
PS fudge factor is equivalent to the inverse of conditional prob- 
ability of finding a region where > Sc of mass scale M[ pro- 
vided that it is totally included inside an isolated region where 
= of mass scale M2. In other words, the PS fudge fac- 
tor is not a constant and depends on smoothing scales M\,M2 
in general. This behavior was also noticed by Nagashima & 
Gouda (1997) using Monte Carlo simulations. Although it has 
been claimed that the EPS formalism provides a good fit to 
the mass function obtained from A^-body simulations for some 
non-Gaussian models (Robinson & Baker 2000), one cannot 
immediately give a justification for the result. In fact, recent 
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numerical simulations of linear density fields showed that the 
EPS formalism does not provide a good fit to the mass func- 
tion for strongly non-Gaussian probability distribution func- 
tions (PDFs) with small variance a < 0.5 which correspond 
to objects with large mass, such as galaxy clusters at present 
(Avehno & Viana 2000). For such small variances, the mass 
function is essentially determined by the abundance of rare den- 
sity peaks which are sensitive to the non-Gaussianity of initial 
fluctuations. It is of crucial importance to understand the role 
of the conditional probability for strongly non-Gaussian density 
fluctuations, especially at large scales. 

In this paper, we study non-Gaussian models in which one- 
point PDF of density fluctuations smoothed on mass scale M 
obeys a distribution with // degrees of freedom, which are 
very simple and widely used as toy models in the literature (e.g. 
Barreiro, Sanz, Martinez-Gonzalez & Silk 1998). The degree 
of non-Gaussianity is characterized by v. A distribution is 
strongly non-Gaussian for small v and converges to a Gaussian 
distribution in the limit v ^ co. 

In §2, we introduce a formalism for computing the mass 
function. In §3, we describe simple models for which a den- 
sity fluctuation at a point obeys a x?- distribution. In §4, we 
explore the property of the conditional probability and com- 
pare the obtained mass functions with those predicted by using 
the EPS formalism. In §5, the effect of mode correlation is dis- 
cussed. In the last section, we sunmiarize our result and draw 
our conclusions. 

2. THEORY OF MASS FUNCTION 

To compute mass functions analytically, PS made following 
assumptions: (1) the overdense region collapses to an virial- 
ized object with mass M when the linear density fluctuation b 
smoothed on mass scale M reaches a critical value be which is a 
function of cosmic time; (2) Each overdense region is indepen- 
dent and described by a spherically symmetric collapse model 
which specifies (5^ (Tomita 1969; Gunn & Gott 1972). Then 
the volume fraction of the collapsing region at initial time with 
mass scale equal to or larger than M is simply given by 



V'(><5oM)= / /?((5;cr(M))c/(5, 



(1) 



where p(d;a{M)) denotes a one-point PDF with variance 
(T^(M) = (Sli) of the initial density perturbation. Now consider 
regions where the density fluctuation S smoothed on mass scale 
Ml exceed 6c. Each region should be totally contained inside 
an isolated collapsed object with mass M2>M\. Then we have 



(2) 



V'(>(5c;Mi)= / ^niM2)PiMi\M2)dM2, 

Jm, p 

where p denotes the mean cosmic density, P(Mi\M2) is the 
conditional probability of finding a region Vi of mass scale 
Ml where i5i = S(Mi) > Sc provided that Vi is totally con- 
tained in an isolated overdense region V2 where S2 = S(M2) = Sc 
(Jedamzik 1995). In this formulation, one can calculate the 
mass function by solving the integral equation (2), once the 
conditional probability P(M\ IM2) and the one-point PDF of the 
smoothed density fluctuations are given. Note that the excur- 
sion set formalism developed by Bond et al. (1991) is essen- 
tiaUy equivalent to this formahsm. Because V2 is an isolated 
region, the conditional probability is given by the probability of 
the first upcrossing of 82 at the threshold 5c when smoothed on 
decreasing mass scales 

P(Mi IM2) = p{5x > 5c\52 = 5c,5iM) < 5c^ M > M2), (3) 



assuming that the spatial correlation in the density fluctuations 
is negligible. Let us call a sequence of density fluctuations 
5(Mi),5(M2),--- in descending order of M as a trajectory of 
a particle. For fluctuations S{M) smoothed by the sharp ^-space 
filter whose phases of Fourier modes (5k are uncorrelated, the 
motion of a particle is described by a Markovian random walk. 
In this case, the conditional probability P(Mi IM2) does not de- 
pend on the state of a particle before the crossing 82 = 5c, 



P{Mi\M2)=p(5r>5c\52 = 5c). 



(4) 



In what follows we assume that equation (4) gives a good 
approximation of the conditional probability for fluctuations 
smoothed by other filters (Yano, Nagashima & Gouda 1996). 
Then we only need to specify a one-point PDF and a two-point 
PDF of the smoothed density fluctuations at a given point. 

If P(Mi IM2) = / does not depend on mass scales Mi and 
M2, then the mass function is described by a formula similar to 
the PS formula, where the PS fudge-factor, /, is related to the 
conditional probability as 



/ = \p(5i >5c\52 = 6c)] ' = const. 



(5) 



Let us first consider the Gaussian models. The bivariate Gaus- 
sian two-point PDF with vanishing means is given by 



PG{5\,52\o\,a2,T)d5id52 = 
1 



277(71(72^1— r-^ 



exp 



2(1 -t2) 



d5id52i&) 



where denotes the variance of 5„ {5}), and r de- 
notes the correlation coefficient, (<5i(52)/(7ifT2. Note 
that pg is scale-invariant, i.e. pg{x\,X2;cf\,(J2,t) = 
a2^ Pc{x\ / (T\ ,X2/ o-2,;l , \ ,t). For fluctuations smoothed by 
the sharp A:-space filter, t = a2/<y\ (see appendix). Then the 
two-point conditional PDF is written in terms of the one-point 

PDF as pciSi \52) = pg(5i , 82)/ PciSi) = Pg(5i - 82; \fo\^-o{). 
Because the Gaussian one-point PDF is also scale-invariant, we 
recover the constant PS fudge factor, f^^ = pci8i > 5c\82 = 8c) = 
(2Tr)-^^^J^exp[-5^/2]d5= 1/2. From equations (2) and (4), 
we obtain an explicit form of mass function. 



niM;5c)dM = - 



fpdi;i>8c;M) 



M 



dM 



dM. 



(7) 



In general, however, the PS fudge factor / is not a constant 
and the mass function n(M) cannot be written expUcitly as in 
equation (7). For instance, for Gaussian fluctuations smoothed 
by other window functions, P(M\ IM2) depends on smoothing 
scales Ml and M2 (Nagashima 2001). 

In the non-Gaussian models, the mass scale dependence 
of f (Ml IM2) must be always taken into account. However, 
Koy ama. Soda <fe Taruya (1999; hereafter KST) claimed that for 
generic non-Gaussian models, the relation p(8\ \ 82) — p(8i — 82) 
holds for Ml <C M2. Consequently, we have 



P{Mi\M2)' 



POO 

■ / Pi5l,(Tl 

Jo 



)dSi 



(8) 



If the one-point PDF is scale-invariant, the above integration 
yields a constant value. Although obtaining a similar rela- 
tion for all scales Mi < M2 is a more complex issue, KST 
and some authors evaluated mass functions by solving equa- 
tion (7) assuming a constant / for various non-Gaussian mod- 
els (Koyama, Soda & Taruya 1999; Robinson & Baker 2000). 
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From now on, we call the PS formalism using the approxima- 
tion described by equations (7) and (8) in evaluating the mass 
function, the extended PS (EPS) formalism. 

Now let us evaluate the validity of the EPS formalism. In the 
limit of vanishing correlation coefficient, t ^ 0, or equivalently 
Ml <C M2, the two-point PDF is written as a direct product of 
one-point PDFs, piS\,52) = pi5\)pi62), which gives 



P(Ml|M2): 



p{Si;(Ti)d6i 



(9) 



On small mass scales with large variance, i.e. ai ^ 6c, the 
lower limit of integration variable 5c can be set to zero as in 
the Gaussian fluctuations smoothed by the sharp fc-space fil- 
ter, leading to a constant /. Hence, the EPS approximation is 
valid for scales Mi <^ M2 and a\ 6c. However, on large mass 
scales, where cri < Sc, such approximation cannot be verified 
except for the Gaussian cases, since the contribution of integra- 
tion of p(6i;ai) in di from to Sc cannot be negligible. There- 
fore, the EPS approximation is not valid for cti < 5c though 
Ml <C M2. This contradicts the validity of the EPS approxi- 
mation for Ml ^ M2 that KST have claimed. In generic non- 
Gaussian PDFs, the relation p(d i\62) ^ p(S i - 62) for Mi <C M2 
does not hold. In fact, the class of PDFs of density fluctuations 
that satisfy p{6i \ 62) ~ pi5i - 62) for Mi -C M2 is very limited' . 
The Gaussian fluctuations smoothed by the sharp A;-space filter 
belong to this very limited class. 

It would be worthwhile to comment on relationship between 
the conditional probability and the excursion set formalism. If 
the conditional probability satisfies pi5i\52) = pi6i-62), then 
the master equation can be reduced to the diffusion equation by 
using the Kramers-Moyal expansion which is used in the ex- 
cursion set formaUsm. However, it is clear that general PDFs 
do not necessarily satisfy the diffusion equation. We need to 
derive a proper two-point PDF that corresponds to the density 
fluctuation distribution function under consideration. 

3. MODELS 

We consider toy models in which the density fluctuation 5m 
smoothed on scale M at a given point obey a xl PDF and we 
assume that the Fourier modes of fluctuations (5k are totally un- 
correlated. The validity of this assumption is discussed in §5. 

Let X; = (xi,-,X2;, • • • ,XMi), i = 1, - ■ ■ ,iy be independent Gaus- 
sian A^-dimensional vector variables. We call the distribution of 
z = (Er=i4'Er=i4'---.Er=i4,) » A^-polnt distribution 
with v degrees of freedom. The one-point xt distribution is 

(i/-2)/2 

p^2 (x; a)dx = 



aT(u/2) 



V X 

la 



exp 



V X 

la 



dx, 
(10) 

where T{x) denotes the F function. 

Next, we derive the two-point PDF p{z\,Z2) of x^ distribu- 
tion with variance (Tj , (T2 and correlation coefficient r. Let us 
consider variables xi,X2 which obey a two-point Gaussian dis- 
tribution pg(xi,X2;si ,S2, e) with vanishing means. By changing 
the variables as zi = x\, Z2 = x^, one obtains 



p^2(zi,Z2;T)dzidz2 ■■ 



27riii2-\/(ziZ2)(l-e^) 



cosh 



X exp 



S2Z1/S1+S1Z2/S2 
2siS2(l-a 



dzidz2, (11) 



where sj = (xf) 
{z\){z2))/<y\a2 = e^. 
is 



= (zi), af = 2sf and r = ((ziZ2) - 
The corresponding characteristic function 



-2iSiS2(slflS2'+*2f2*l') 



■(?-)ht2. 



-1/2 



(/)^2(fi,f2) = [l- 

(12) 

The characteristic function for 77^2 (zi,Z2;o'i,cr2,T) is given by 

0^2 = ii'xi^" ^^^'"^ (^') ~ V^fei) = '^1 = and t = e^, 
since each variable Zi is written as a sum of independent ran- 
dom variables. From two-dimensional Fourier transform of 
(px2(ti,t2), we finally obtain the two-point PDF of xl distri- 
bution. 



Pxlizi,Z2',cri,cr2,T) 

Z1Z2 \ " 



1 



exp 



rj'zi+0-2'z2 



xL 



for 0<r<l, 



= < 



-Z1Z2 



Qa2r(f) 



exp 



'•l'zi+<T2'z2 



(13) 



Xjji_i 



TZIZ2 



aa 



Z1Z2 



for - 1 < r < 0, 



^2(r(^))2 



exp 



^l'zi + (72^22 



for T = 0, 



where a = I -t, a = ^jT/Jxa^Jv and y„ and /„ are the Bessel 
and modified Bessel function of the first kind, respectively. The 
one-point and the two-point PDFs are scale-invariant and 
extend from to 00 for each variable. Because we assume that 
the PDFs of a density fluctuation have a vanishing mean, we 
wiU use off-centered PDFs, p*^^ (z) = P^l^z-^ (z))- P^' (^i'^^) = 
p^2 (zi + (zi ) , Z2 + (z2)) which are also scale-invariant in the fol- 
lowing analysis. 

It is cumbersome to calculate the three-point PDF explicitly. 
For simplicity, we only consider the case of = 1 (see also 
Sheth 1995). In terms of the covariance matrix for the origi- 
nal trivariate Gaussian PDF, 



cov = 




(14) 



the explicit form of the trivariate Xi PDF can be written as 

/'x?fel'22,Z3;0'l,f2,0'3,(Ti2,a-23,fl3) 



1 



:exp 



C11Z1+C22Z2+C33Z3 



2(27r)3/2^det(cov)ziZ2Z3 

X (^exp(-Cl2 V2122) COSh(C23 V^iis + C13 -y/ziZs) 
+ exp(Ci2v^ZlZ2)cOsh(-C23\/Z2i3+Cl3\/ZlZ3)), 

where cry = 2i?- for i ^ j, af = ls% and = (cov"')y. 



(15) 



' Although KST claimed that for Mi <C M2, the relation p{Si \S2) ~ p(Si -S2) holds for generic non-Gaussian models, the statement is incorrect, 
correct if the cvmiulants satisfy a relation {5"^5'^"')c = {62)': ^^'^ sXln>m rather than moments. For Gaussian fluctuations, all the cumulants {5"5'^ 



for n -- 



The statement is 
')c vanish except 

2 (the means are assumed to be zero). However, for generic non-Gaussian fluctuations, the cumulants do not necessarily vanish. A condition for moments that 



has been used in KST as an ingredient of the proof {S'^S^ "*) = (Sj) is also incorrect for n>2. For instance, {S\5\) - 



2(t\< 



-\ for T = 



0. 
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It should be noted that the models considered in our analy- 
sis are not exactly identical to the "xt field models" (Peebles 
1999; Scoccimarro 2000) in which the initial fluctuation itself 
is described by a xt field. In contrast to the Gaussian mod- 
els, the PDF of a smoothed x^ held deviates from the original 
PDF in general (Avehno & Viana 2000). If one would hke to 
explore more realistic models, one should take into account the 
dependence of the PDF on smoothing scale. Fortunately, for xt 
PDFs with p larger than 10 smoothed by the Gaussian filter, it 
is known that the departure from the original PDF is not signifi- 
cant although the departure is noticeable in the case of the sharp 
A:-space filter (Avehno & Viana 2000). Hence, we can expect a 
similar result for some particular choice of smoothing filter. 

4. SOLVING THE CLOUD-IN-CLOUD PROBLEM 

In order to study the characteristics of x^ models on mass 
function, we first estimate the conditional probability 

P(Mi|M2)= / p*^2iSi\6c;au(T2,T)d6i (16) 

J 6c 

for various kinds of smoothing filters, such as the sharp A:-space, 
the Gaussian and the top-hat filters. 

First of all, we consider fluctuations smoothed by the sharp k- 
space filter Wr(A:) = 6B?0{-k /R—k)/-!:. The relation between the 
mass M and the smoothing scale R is given by M = dR^pj'K. The 
correlation coefficient is t = ajjox = {M\/M2f"^^^^^ (see ap- 
pendix) provided that the power spectrum of the initial density 
fluctuations has a form V{k) oc k", where n denotes the spectral 
index and n > -3. The critical value 5c is 1.69 independent of 
mass scale for the spherical collapse in the Einstein-de Sitter 
universe. For simplicity, we assume 5c = 1.69 in the following 
analysis. As we have argued in §2, for Mi <c M2 and <j\ 5c, 
the EPS approximation gives a correct value, 
1*00 

r\EVS) = j p*^2j5uai)d5i 

_ T(iy/2,u/2) 
T{iy/2) 

= PiMi\M2) forM2>Mi and(7i >^c, (17) 
where r(x,y) denotes the incomplete Gamma function. For 
other parameter regions, however, the EPS approximation is 
not always correct, especially for strongly non-Gaussian PDFs. 
From Figure 1, one can see that the inverse of the EPS factor, 
, deviates from the correct conditional probabiUty P(Mi IM2) 
for a region Mi <C M2 and <j\ < 5^, or equivalently Mi > M*, 
where cr(M*) = 1. For fluctuations smoothed by the sharp k- 
space filter, the trajectories are described by a Markovian ran- 
dom walk. Therefore, the chance of upcrossing at the criti- 
cal value 5c is almost equivalent to the chance of downcrossing 
at 5c, i.e. P(Mi\M2) = 1/2. Therefore, in the neighborhood 
of diagonal fine Mi = M2, we have P(Mi IM2) ^ 1 /2. On the 
other hand, in the Gaussian limit ly — > 00, we have /"'(EPS) = 
T{v/2,v/2)/T{y/2) 1/2. Consequently, for weakly non- 
Gaussian PDFs smoothed by the sharp A:-space filter such that 
/"'(EPS) 1/2, it is natural to expect P{Mi\M2) - 1/2 on 
all scales Mi < M2. Thus, for weakly non-Gaussian PDFs 
smoothed by the sharp A:-space filter, the EPS formalism gives 
a good approximation of /. 

Next, we consider the property of fluctuations smoothed by 
other filters for which the trajectories of density fluctuations 
cannot be described by a Markovian random walk because of 
the Fourier-mode correlation. In this case, equation (4) does 



not give an exact result. Let us consider a trajectory that firstly 
crosses 5c upwards at M2 and ends at i5i > 5c and Mi. For 
M2 Ml, the trajectory can be well approximated by that of a 
Markovian random walk, since the correlation length is almost 
negligible compared with the length of the whole trajectory. On 
the other hand, for M2 ~ Mi, the trajectory can be well approx- 
imated by a monotonically increasing function in decreasing 
value of M, or increasing value of a, since the correlation coef- 
ficient is almost equal to unity. Then the conditional probabihty 
is approximately given by 

P{Mi\M2)^ j p(5i\52 = 5c,d52/dM2<0)d5i. (18) 

At large scales, 5i 5c, or equivalently Mi S> M,, the prob- 
ability of upcrossing at 5c is very low. In other words, the 
probabihty of 052/ dM2 > at 52 = 5c is almost zero. There- 
fore, the extra condition 852/ dM2 > is not necessary at large 
scales. As shown in Figure 2, when the Gaussian filter is used, 
f (Ml IM2) increases compared with that for the sharp fc-space 
filter owing to the Fourier-mode correlation. At large scales 
Ml M*, f (Ml IM2) ~ 1. Consequently, the chance of down- 
crossing at 5c is almost zero. Thus it is reasonable to conclude 
that equation (4) still provides a good approximation of mass 
functions, at least, at large scales even for other filters. On the 
other hand, in the region Mi < M* and Mi < M2, P(Mi IM2) 
is close to the value in the case of the sharp A:-space filter. This 
suggests that the approximation of the Markovian random walk 
is good in that parameter region. 

Substituting the variances ai, 02, and the correlation coef- 
ficient r for each type of smoothing filter (see appendix for 
derivation) into the conditional probabihty equation (4), one 
can compute the mass function by solving the integral equa- 
tion (2). In Figure 3, we show the multiplicity functions F(M) 
for the various filters and for the EPS formalism in the case of 
Xi PDF with one degree of freedom. The multiplicity func- 
tions derived here decrease at large scales compared with those 
of the EPS prediction. As we have argued in §2, the actual 
value of P{M\ IM2) = p{5\ > 6c\52 = 5c) is larger than the value 
/"' predicted by using the EPS approximation at large scales, 
which explains a decrease in the multiphcity function F{M). To 
compensate the deficit at large scales, F{M) increases at small 
scales. Regarding the dependence of F{M) on the spectral in- 
dex n, one can see in Figure 3 that a relative increase of F{M) 
at small scales is much prominent for larger n. This is because 
F(M) decreases rapidly at smaller scales for a larger value of n, 
or for a much blue spectrum. Similar behavior is also observed 
for the Gaussian models (Nagashima 2001). 

In Figure 4, we show, for PDFs with various degree of 
freedom, dgia) /da which is defined as gia) = g{M{(j)) for 
which g(M) = F(M')dM' is the fraction of collapsed objects 
above a smoothing scale M. At large scales, a <l, dg(a)/da 
is significantly increased compared with those corresponding to 
the Gaussian PS mass function. This is because the xt PDFs 
have a broad tail toward large 5, especially for strongly non- 
Gaussian PDFs with small value of Even in the case of 
v = 50, one can still observe a clear difference from the value 
corresponding to the PS mass function for fluctuations with 
small variance < 1 . In other words, the amount of rare density 
peaks at large scales is very sensitive to the non-Gaussianity 
of the initial fluctuations. In the Gaussian limit, ly ^ 00, the 
EPS formalism correctly reproduces dg{a)/da for fluctuations 
smoothed by the sharp fe- space filter. As the degree of freedom 
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V decreases, a deviation from the EPS prediction becomes no- 
ticeable. It is clear that the EPS formalism overestimates the 
number density of dark halos on large scales cr < 1, especially 
for strongly non-Gaussian PDFs. For fluctuations smoothed by 
the Gaussian and the top-hat filters, such a deviation is promi- 
nent even in the case of a weakly non-Gaussian PDF (y = 50) 
owing to the correlation between Fourier-modes. Similar re- 
sult has been obtained by AveUno & Viana (2000) using Monte 
Carlo simulations for smoothed xt fields. 

Note that our analytic calculation of dg{a)/da for = 10 
and IT < 1 agrees well with the previous numerical results by 
AveUno & Viana (2000) when the Gaussian filter is used. The 
result is natural, since the departure from the original one-point 
PDF is not significant in this case. However, for fluctuations 
smoothed by the Gaussian filter with sUghtly large variance 
l<cr<5, our analytic values exceed their numerical values. We 
might be able to recover their numerical results if we use the im- 
proved approximation described by equation (18) of the condi- 
tional probability instead. Nevertheless, for the purpose of con- 
straining non-Gaussianity of initial fluctuations using the abun- 
dance of high-redshift clusters (Matarrese, Verde & Jimenez 
2000), our approximation of the conditional probabihty may be 
sufficient in practice, as mentioned. 

5. EFFECT OF MODE CORRELATION 

So far, we have assumed that different Fourier-modes of x?. 
fluctuations are uncorrelated. However, if we consider field 
models defined as 0(x) = V"? (x) where V'i(x) are indepen- 
dent homogeneous and isotropic random Gaussian fields whose 
Fourier-modes are uncorrelated, Fourier-modes of a x^ field 
are no longer uncorrelated. The signature of correlation ap- 
pears as non-vanishing skewness , kurtosis, or higher-order cu- 
mulants although the two-point correlation has the same form 
as the Gaussian one (Scoccimarro 2000). In this case, a hypo- 
thetical particle with position 5{M) cannot be described by the 
Markovian random walk even if we smooth the fluctuations by 
using the sharp-k space filter Therefore, it is of crucial impor- 
tance to estimate the effect of correlation. 

In order to check the validity of our calculation, we com- 
pare the conditional probability Pt, = p(5\ > 8c\52 = S^Sj, < 
Si), Ml < M2 < M3 specified by a three-point and two-point 
PDFs to P2 = p(Si > Sc\52 = Sc) specified by a two-point and 
one-point PDFs. In fact, the actual condition to be satisfied is 
equation (3), which is equivalent to P„,n — > 00. P3 is expected 
to be more accurate than P2 for estimating the conditional prob- 
ability P(Mi IM2). 

In what follows, we only consider v = I case, for which 
the effect of mode correlation owing to the non-Gaussianity 
is most significant. In the limit 0-3 0, the correlation co- 
efficients vanish ^ 0-3/0-1 0, which gives p(Si\S2,Sj,) = 
p(5\,52)p(h)/p(52)p(h) = p{5\\52)- Thus, if M3 is large 
enough compared with M2 and M\, the effect of mode corre- 
lation is negligible. On the other hand, for M3 ~ M2, or equiv- 
alently, (73 ~ 02, the probability of having a trajectory 5\ > Sc 
becomes large owing to the correlation in the fluctuations on 
different scales. Note that in the hmit 03 0^2, P3 is equivalent 
to the right-hand-side of equation (18) assuming that (52(o'2) is 
differentiable at 02- From equation (15), one can easily calcu- 
late flie ratio P3/P2- 

First of all, we consider the case of the Gaussian filtering. It 
turns out that the ratio /^3 7/^2 is < 1.5 for ai = 2, and < 1 .2 for 
(Ti = 0.5 if the spectral index is n = -2 (see figure 5). At large 



scales (J < 1, a decrease in the multiplicity function F(M) is 
expected to be less than ~ 30 percent. 

In the case of the sharp-k space filtering, we have the scale- 
invariant relations in the covariances: an = o\,a23 = 0-^,0-13 = 



C7|. Then we obtain 

P^JiZl\Z2,Z3) 
= pg(Zl|Z2) 



4) 



= exp 



Z1+Z2 



cosh 




(19) 



, (20) 



where of = 2sj, which gives P3 = P2 irrespective of the value 
of the spectral index n as in the case of the Gaussian fluctua- 
tions. Thus, the effect of three-point correlation is completely 
negligible if the fluctuations are smoothed by using the sharp-k 
space filter. This relation might also hold for other PDFs 
with y > \, since the effect of correlation is less significant. 
Presumably, all the scale-invariant PDFs have this property. 

To summarize, the effect of mode correlation is less signif- 
icant for fluctuations at large scales o < 1 irrespective of the 
type of filtering. Inclusion of an additional condition ^3 < dc 
improves the reUance of calculation for strongly non-Gaussian 
fluctuations at smaU scales a > I. 

6. CONCLUSIONS 

In this paper, we have demonstrated that the EPS approach 
has a drawback in describing the number of collapsed ob- 
jects, especially for strongly non-Gaussian density fluctuations. 
Based on the formalism developed by Jedamzik, and Yano, 
Nagashima and Gouda properly taking into account the scale 
dependence of correlation between objects at different scales, 
we analytically calculated the mass function for various 
models and found a deviation from those predicted by using 
the EPS formalism, especially noticeable for strongly non- 
Gaussian models: a decrease at large scales and an increase 
at small scales in the value of multiplicity function. The re- 
sult may affect some recent studies of constraints on non- 
Gaussian models using the cluster abundance at different red 
shifts and the correlation length of galaxy clusters (Chin, Os- 
triker & Strauss 1998; Koyama, Soda & Taruya 1999; Robin- 
son, Gawiser & Silk 2000). Our results are similar to those by 
Avelino and Viana (2000) based on Monte-Carlo simulations of 
non-Gaussian xt fields. At intermediate scales, 1 < cr < 5, the 
deviation from the EPS prediction is not prominent. It seems 
that the result is consistent with those from A^-body simulations 
for various non-Gaussian fields (Robinson & Baker 2000). It 
would be interesting if larger A'-body simulation could be car- 
ried out and find out a deviation of mass function from the EPS 
prediction for objects with very large mass (low-cr) or for those 
with very small mass (high-cr). 

In order to vindicate that the generalized PS formalism works 
in estimating the mass function for non-Gaussian models, we 
should take various kinds of effects into consideration: effects 
of non-spherical collapse (for Gaussian fields, see Sheth et al 
2001), ambiguity in mass-smoothing scale relation (Bardeen 
et al 1986, Peacock & Heavens 1991), and conditions of ob- 
jects surrounding by an isolated dark halo. The last issue is 
relevant to the cloud-in-cloud problem. Although we have dis- 
cussed about a prescription for incorporating the condition of 
upcrossing at the critical value 6c, we did not explicitly con- 
sider the effect of spatial correlation owing to the finite size of 
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halos in evaluating the conditional probability P(Mi\M2)- For 
Gaussian models it is known that the effect of spatial correlation 
almost cancels out the filtering effect, recovering the original 
PS mass function, particularly in the case of the top-hat filter 
(Nagashima 2001). It is of very importance to check whether 
such a cancellation occurs for non-Gaussian models. The other 



issues left untouched should be addressed in our future work. 
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cussion and comments. We would also like to thank the anony- 
mous referee for useful suggestion and comments. 



APPENDIX 
CORRELATION COEFFICIENTS 

In this section, we derive the formulae of correlation coefficients t for density fluctuations smoothed by three types of filter, 
namely, the sharp A:-space filter, the Gaussian filter and the top-hat filter. 

Let us consider a mass density fluctuation (contrast) S(x) = ^^y^, where p(x) denotes the mass density at a point x and p is the 
mean cosmic mass density. Using a window function WsCr = \x\), (5(x) can be smoothed on scale R 

SRix) = JwRi\x-x'\)Six')dx, (Al) 

and in Fourier space, 

= TT^ / WRimk)e''^-^dk, (A2) 

(27r)J J 

where St and WR(k) are the Fourier transforms of 5^ and W^Cr), respectively. Here we choose a normalization of the window function 
asW^R(0)= 1. 

If the fluctuations are homogeneous, then the two-point correlation is written in terms of the power spectrum V{k) as 

(^k, = (27r)'fe(ki -k2)VikO. (A3) 

which gives 

((&)') = 7Z-T3 / 47rieW^ik)Vik)dk. (A4) 

The correlation coefficient is given by r = {Sr^Sr^) / cfr^ cfr^ where 

1 f°° 

{SrM) = 7T^ 4TTieWR,ik)WR,ik)Vik)dk. (A5) 

Note that aU the equations in this appendix except for (A3) hold for fluctuations whose Fourier transforms are correlated. The mass 
of objects smoothed on scale R can be defined as 

M= j pWRix)dx. (A6) 
The correlation coefficients for the three types of filters are written as foUows. We assume that Vik) oc fc". 

1. Sharp A:-space filter 

3 

WRir)= -pr^ (&m kcr-kcr cos kcr), (A7) 

ky 

WRik)=^9ikc-k), (A8) 

where kc is the cut-off wave number, kc — R^^, and 9(x) is the Heaviside step function. Note that there is an ambiguity in the 
relation between kc and R. Here we define it as kc = n/R that gives M = 6R^p/tt which has been widely used in the literature. 
If we choose kc = (97r/2)^/^//?, then we have M = AnR^p/S. For the former definition, the variance of density fluctuation is 

and the correlation coefficient for Mi < M2 is r = auJcTM, = (M2/Mi)-<"+3)/6 = (/?2//?i)-(«+3)/2 

2. Gaussian filter 

WRir)=expi-r^/i2R^)), (AlO) 

WR(k) = {l-Kfl^R^ exp(-k^R^/2). (Al 1) 

The mass of objects smoothed on scale RisM= (2tt)^^^R^p and for n > -3, 

a^oc(2^)-2/?-3-"r((3 + n)/2), (A12) 

9(n+3)/2/^-4/(3+n) -4/(3+«)._(3+„)/2 

r = = ^ '- . (A13) 
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3. Top-hat filter 



WR{r) = e{\-r/R), 



WR(k) = AnK 



(sin A:/? cos kR 



The mass of objects smoothed on scale /? is M = AwR^p/S. If n is not an integer, for -3 < n < 1, 

(tIi cx ^-^(2/?r(n-2)(n+l)sin(«7r/2)r(-3+n), 



(A14) 
(A15) 

(A16) 



r=^^S^( [iR^+R2f-''-(R2-Rif-" m-3 + n) + T(-2+n)) 



+ (^{Ri +/?2)'""-(/?2-/?i)'"")/?i/?2r(-l +n)j a^'af. 



For integers -3 <n < 1, 



and 



r = < 



/ TT 

6 ^" = «^ 



(2/?)-" X < 



I (« = -l) 



^ 60 (" = -2) 




, (n = 0) 



2RiR2{R{ +/?2) + (^2 -^i) In 



/?l+^2 



5Rl-R\ 



4^3 -v/ra=^(5at-c72^), (n = -2). 



(A17) 



(A18) 



(n = -l) 



(A19) 
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Fig. A1. — Plots of the conditional probability p(Si > Sc\S2 = Sc;cr(Mi),a(M2),T(Mi,M2)) for xt PDFs smoothed by the sharp i-space filter, n denotes the 
spectral index. The accompanying palettes show the relation between the gray level and the value. M, is defined as it(M,) = 1. We set Sc = 1.69. 
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Fig. A2. — Same as Figure 1, but for the Gaussian filter. Note that the contour level is from to 1. 




-1-2 2 



INOUE and NAGASfflMA 



9 




-6-5-2 ?. 4 

Log[M/M*] 




-e -4 ~2 

Log[M/M*] 



Fig. A3. — Plots of multiplicity functions F(M) =Mn(M)/p for fluctuations obeying the Xi PDF smoothed by various filters. (a)n = -2. (b)n = 0. The solid, the 
dashed, and the dashed-dotted lines correspond to F(M) for fluctuations smoothed by the sharp fc-space filter, the top-hat filter, and the Gaussian filter, respectively. 
The dashed-double-dotted line represents the value predicted by using the EPS formaUsm. M* is defined as it(M«) = 1. In the case of n = 0, the multiplicity function 
F{M) for fluctuations smoothed by the top-hat filer is exacUy the same as that smoothed by the sharp k-space filter. 
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Fig. A4. — Plots of dg/da for fluctuations that obey three kinds of PDFs (u = 1, 10 and 50). The spectral index is assumed to be n = -2. The solid, the 
dashed, and the dashed-dotted hnes correspond to F(M) for fluctuations smoothed by the sharp /:-space filter, the top-hat filter and the Gaussian filter, respectively. 
The dashed-double-dotted line represents the value predicted by the EPS formalism. The dotted line corresponds to the value predicted by the PS formalism for the 
Gaussian PDF. In the case of the sharp-k space filter, dg/da does not depend on n. 
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Fig. A5. — Effect of mode correlation in tiie case of Gaussian filtering for Xi fluctuations with the spectral index n = —2. A deviation from 1 in the value of P3/P2 
gives a relative error in for estimating the conditional probabflity P(Mi IM2) assuming that the effects of four-point or higher-order correlations are negligible. 



